Hubbard ring: currents induced by change of magnetic flux 
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We investigate currents in a quantum ring threaded by a magnetic flux which can be varied in an 
arbitrary way from an initial value <j>i at time U to a final value <f>f at time tf. Dynamics of electrons 
in the ring is described by the Hubbard and the extended Hubbard models. We demonstrate that 
time dependence of the induced current bears information on electron correlations. In the case of the 
Hubbard model with infinite on-site repulsion we prove that the current for t > tf is independent of 
the flux variation before t / . Additionally, this current is fully determined by a solution of the initial 
equilibrium problem and the value of <f>f. Apart from mesoscopic rings our results pose important 
implications for designing of quantum motors built out as the ring-shaped optical lattice. 
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Time-dependent manipulation of quantum states in 
nanosystems is an important problem directly related 
to future applications both in the context of quantum 
control [J] and the reduction of decoherence Q- There 
is a natural ground for implementing such ideas: meso- 
scopics, where quantum effects play a crucial role Q. 
Unfortunately, analysis of quantum systems affected by 
external time-dependent force or field F(t) is extremely 
difficult and only very few models are exactly solvable. 
The best known examples concern a driven quantum os- 
cillator and a two-level system in circularly polar- 
ized magnetic field 0, 0] • In many cases the solvability is 
related to certain symmetries of the system. Moreover, 
there are interesting regimes of strong external driving 
when the linear response theory cannot be applied. For 
periodically driven systems one may apply the Floquct 
theory however, in many cases this general approach 
provides only approximate numerical results. Another 
problem concerns the possibility of experimental verifica- 
tion of theoretical predictions. Here a pertinent question 
arises: how does the final (generally non-equilibrium) 
state depend on a particular shape of -F(t)? For the pur- 
pose of practical applications it is desirable to find sys- 
tems which are robust against small temporal changes of 
F(t), e. g., originating from imperfect realization of the 
assumed conditions. A realistic example will be pointed 
out in this work. 

We consider one-dimensional (ID) mesoscopic systems 
of ring topology threaded by a magnetic flux <fi. Dynam- 
ics of electrons moving in the ring is described by the 
Hubbard and the extended Hubbard models. We pro- 
pose a simple scheme involving variation of the magnetic 
flux to manipulate currents: it is changed from its ini- 
tial value <fi(t) = <f>i for time t < t$ to the final value 
4>(t) = <fif for t > tf. We show that for vanishingly small 
and infinitely strong many-body interactions the result- 
ing current does not depend on the way the magnetic flux 
is modified or switched on. In other words, the current 
for t > tf is entirely determined by the solution of the 
equilibrium problem at t = tj. Although the method of 
reasoning we apply for non-interacting fermions is rather 



trivial, it nicely illustrates the general method that is ap- 
plicable also to a non-trivial case of correlated electron 
systems. Our analysis can be applied either to nonsuper- 
conducting mesoscopic rings or to rings built in the op- 
tical lattice setup [10(. The difference in energy scales in 
both systems shows up mainly in different time scales of 
the external driving. 

We start with the Hamiltonian of non-interacting par- 
ticles in the ring threaded by a magnetic flux. It is a sum 
of one-particle Hamiltonians 



H{t) = 7T 
2m 



(1) 



where L is the circumference of the ring and e is a charge 
of the particle. The current operator is related to the 
momentum observable in the following way 



l(t) = = -L- 

1 ' dcf>(t) ml 



(2) 



Now, let us assume that the magnetic flux cj) is varied in 
an arbitrary way from an initial value <fn at time ti to a 
final value 0/ at time tf. One can explicitly extract the 
time dependent part of the current operator 



I(t)=Ii+AI(t), 



AI(t) = 
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The averaged current flowing in the ring is determined 
by the relation 



(/(*)) =Tr [/>(*)/(*)] 



(6) 



where p(t) is a density matrix of the system. Its time 
evolution is determined by the von Neumann equation 
ihp(t) = [H{t),p{t)]. Note that from Eq. ||T]) it follows 
that 



[H(t 1 ),H(t 2 )} = 



(7) 
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for arbitrary t\ and t^- In consequence, a solution of the 
von Neumann equation has the form 



p(t) = exp 



- / H(s)ds 



to 



p(i )exp 



- / H(s)ds 



.(8) 



If p(to) commutes with the Hamiltonian H{t) for any 
t > to then p(t) — pita). The latter requirement is not 
very restrictive. It is fulfilled by the Gibbs state p(to) oc 
exp [— PH(to)] and, from the experimental point of view, 
seems to be the easies and the most natural choice of the 
initial preparation. Let us take t$ = ti and pi = p{ti). 
Then the averaged current reads 



</(*)> = Tr [p(t)I(t)} = Tx\pili] + AI(t), 



(9) 



where the first term on r.h.s. is the initial equilibrium 
persistent current and the second term is the current in- 
duced by the time-dependent component AI(t). As the 
latter quantity is independent of the initial state, one 
can easily calculate the current induced in a system of n 
non-interacting particles AI n (t) = nAI(t). It is instruc- 
tive to compare AI n (t) with the amplitude of persistent 
currents Ipc at zero temperature . One finds 



<AJ„(i)> MU)-<t>{t) 



= 2 



L PC 



(10) 



where Ipc — vpe/L, vp = tmnj (ml) is the Fermi veloc- 
ity and 0o = h/e is the flux quantum. 

Let us notice remarkable properties resulting from 
Eq. tfTO")) : (1) The averaged current (I(t)) depends on the 
flux only at the same instant of time t, i.e., the current is 
independent of the way how the magnetic flux is switched 
on; (2) (I(t)) is fully determined by the solution of the 
initial equilibrium problem, i.e. by the equilibrium mean 
value Tr[/9jij]; (3) one can induce currents which are sig- 
nificantly larger in amplitude than the persistent currents 
provided that <j)(t) — 4>{ti) 3> fa- 
It is known that electrons in ID systems are almost 
always strongly correlated [12j. In that sense, the free 
electron approximation applied to ID quantum rings is 
disputable. In the following we discuss to what extent 
the conclusions derived for the model of noninteracting 
particles are applicable to more realistic systems of cor- 
related particles. Here we consider a ID ring described 
by the 'standard model' of correlated electrons, i.e. by 
the Hubbard model [IlEl, 



H H {t) 



-J 



3,° 



a j+l,a a 3,cr 



h.c 



We choose J as the energy unit, whereas time and current 
will be expressed in units of r = H/ J and Iq = 2nJ /Nfa, 
respectively. 

Before we carry out discussion based on analytical re- 
sults, it is instructive to inspect numerical studies. Fig. 
1 illustrates the time-dependence of the average current 
(I{t)) obtained from Eq. © by a solution of the von 
Neumann equation for p(t) under an initial condition 
being the equilibrium state pi oc exp [— /3Hn(ti)} with 
(3 — > oo. The presented results have been obtained for 
a ring consisting of N — 6 sites with various numbers 
of spin up n-f and spin down n± particles. Apart from 
the case n-\ + = N (when the system is insulating) 
the qualitative results are independent of the system 
size and the number of particles. We find that in the 
small and large U / J-limits, a dc current is observed for 
t > tf. Its amplitude is independent of the rate of flux 
variation and is greater than the equilibrium persistent 
current. For moderate values of U/J, the current dis- 
plays time-oscillations. However, its average over time 
is non-zero and the dc component can be detected. Fre- 
quency of oscillations depends on electron correlations: 



<Kt)>/i 

or 



<l(t)>/l 



U=J 





H'Slf 



V'H'I 



!' 1 



W»»+HMM»» ; - 'Ht-tttl 2 

t/T 100 50 t/T 100 



(11) 



where J is the hopping integral, U is the on-site Coulomb 

|, j. For the ring con- 

2ir(f>(t)/N(f)Q. The current oper- 



repulsion, 

sisting of N sites <f>(t) 
ator reads 

2vrJ 

3,v 



a j,a a 3,a anC l a 



I(t) 



h.c. 



(12) 



FIG. 1: (color online) The averaged current {I(t)) induced 
by change of the magnetic flux <f>(t) from cj>(ti — 0) = to 
4>(tf) — 4>o with the change rate <f>(t) = const for t £ (ti,tf). 
The results are obtained for iV = 6, rif = 2, ny — 2 (left 
panels) and N = 6, n-f = 3, rij. = 2 (right panels). Continuous 
(red) lines show results for tf — r and dashed (blue) lines for 
tf = 25t. The horizontal dotted lines show maximal values 
of the equilibrium persistent current. 
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for stronger correlations, i.e., for larger U/J, frequency 
of the current is higher. On the contrary, its amplitude 
decreases as U/J increases. The amplitude of oscilla- 
tions is more sensitive to the rate of the flux changing 
= OAf — — slow changes of the flux result 

in small amplitudes of the current oscillations and vice 
versa. The detailed analysis of currents (e.g., regular vs. 
chaotic behavior) in this intermediate regime will be pre- 
sented elsewhere. 

The numerical results suggest that the conclusions for- 
mulated for free particles [(i)-(i?)] hold true also for the 
system described by the Hubbard Hamiltonian (fTTj) with 
U = or U/J = oo. The former case (U = 0) is again 
trivial, since in the Bloch representation one gets 

H H {t) =-J2 2 Jcos(fc - 4>{t))al a a k ^. (13) 

fc,<7 

As Eq. ([7]) holds true for the above Hamiltonian, p{t) = 
p(ti) = const and (I(t)) is independent of the magnetic 
flux <j)(t') for t' < t. In the following we prove that it also 
holds true for U/J — oo. In the case of infinitely strong 
Coulomb repulsion one can rewrite the Hamiltonian pip 
in the form 

H H {t) = H L (t)+H R (t), (14) 
H L (t) = H R (t) = -Jc^Pj2a] +1 ^P, (15) 

where the operator P = Il|L 1 (l — n j,T n j,l) projects 
out states with doubly occupied sites. It is clear that 
[H L (h),H L (t 2 )} = [H R (tx),H R {t2)] = 0. Then, in order 
to prove that Eq. ([7]) holds true, it is enough to show that 
[H L (ti),H R (t 2 )]=0. One finds 

[H L (t 1 ),H R (h)] = J 2 e<^-^)(A-B), (16) 
A = Pa \+i,a a i,TPa) ttL a j+ i tli P, (17) 

i,j,a,)Ji. 
i,j,a,/i 

We have introduced the notation P = P to distinguish 
different positions of these operators. Analyzing the op- 
erator A one can note that because of the presence of 
the projection operators P one can neglect P unless the 
hopping a\ +1 a ai^ removes double occupancy generated 

by at aj+x,^. Similar method of reasoning applies to the 
operator B. Therefore, in Eqs. (fTT)) and (TT5)) one can re- 
place P with (1— Sij)+8ijP. Note that a^-aPa} a \ip) =0 
for arbitrary state Taking into account these prop- 
erties one gets 

A-B = (1 - Sij) p Wi + i,<r*i,°i »L°i+ 1 ./*] p 

+ Y Pa i+l,a a h<T Pa i,a a i+h<rP 
ia 

-^2 Pa i,a a i+h<rPal +l!a ai i(T P (19) 

ia 



The first term vanishes because the commutator is pro- 
portional to Sij. Now, the projection operators P in the 
second and third terms can be replaced by 1 — rii- a and 
1 — ni+x,- CT j respectively. Then, one gets 

A — B = Pj^ni+iA 1 ~ - 
-Pj^niA 1 ~ m+i,<r)0- - n i+ i,-a)P- (20) 

ia 

One can see that A — B is expressed as a difference of two 
operators (first and second lines in the above equation). 
The first one counts how many occupied sites succeed 
empty sites, whereas the latter one counts how many oc- 
cupied sites precede empty sites. In a ring-shape system 
these numbers are equal for an arbitrary state. Hence 
A-B =0, next [H L (ti), H R (t 2 )] = and then Eq. (J7J 
holds true. 

It is important to compare our method with the stan- 
dard approach to ID Hubbard model via the Bethe 
Ansatz [14J. The latter one provides a solution for all 
interaction strengths and band-fillings and we refer to 
Ref. [l5| for a comprehensive review on the equilib- 
rium properties of the ID Hubbard model. In partic- 
ular, the elementary excitations are expressed in terms 
of holons and spinons which, in general interact and are 
not independent (Ts| . Then, analysis of correlation func- 
tions within the Bethe Ansatz is by far not straightfor- 
ward. The complete charge-spin separation over all en- 
ergy scales occurs only in the case U — » oo what signifi- 
cantly simplifies calculations of the correlation functions 
[l6T| , Although our approach goes beyond the the linear 
response theory and does not relay on the Bethe Ansatz, 
the results for U — > oo can be interpreted in terms of 
the charge-spin separation. As the vector potential cou- 
ples to charged orbital degrees of freedom, in the case 
of the full spin-charge separations the system responds 
to time-dependent flux in the very same way as a sys- 
tem of non-interacting fermions. Within this conjectural 
interpretation one comes to a conclusion that observa- 
tion of currents induced by time-dependent fluxes gives 
important information on the spin-charge separation. 

Finally, we verify whether the above discussion can 
be generalized for a system that cannot be solved via 
the Bethe Ansatz. For this sake we have carried out 
numerical calculations for the extended Hubbard model 
with the Hamiltonian [ItJ 

n ha n j+1 ^. (21) 

3 a A* 

Fig. [2] shows similar results to those presented in the 
left column of Fig. [1] but calculated for the extended 
Hubbard model with V = OAU. We are unable to prove 
that the general properties of (I(t)) derived for the non- 
interacting system hold true also for the extended Hub- 
bard model in the limit of strong interactions. However, 
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FIG. 2: The same as in the left column of Fig. [T]but for the 
extended Hubbard model with V = OAU. 



case of the ac-driven quantum motors (considered, e.g., 
in Ref. [l(|) the systems with U ~ J should provide the 
best performance. Generation of significant currents in 
a system with U / J <C 1 or U/J ^> 1 is possible provided 
</>(£/) — <j){ti) ~ 0o- In the latter case, magnitude of the 
current is independent of the way how the magnetic flux 
is modified. This feature may facilitate the experimen- 
tal realization. Summarizing, our proposal expose new 
prospects of inducing, controlling and manipulating of 
currents in nonsuperconducting quantum small systems 
of ring topology. 

The work supported by the Polish Ministry of Sci- 
ence and Higher Education under the grant N 202 131 
32/3786. 



a comparison of Figs. Q] and [2] strongly suggests that it is 
actually the case. Fluctuations of {I it)) for t > tf grad- 
ually extinguish when the interaction becomes stronger 
and the magnitude of these oscillations decreases when t f 
increases. Simultaneously, for U 3> J the magnitude of 
the current for t > tf becomes almost independent oft/. 
We have carried out similar calculations for V = J (V is 
independent of U) and found that the above conclusions 
concerning the role of U remain valid. 

In conclusion, we have analyzed currents induced by 
temporal changes of an external magnetic flux piercing a 
quantum ring. On the one hand, experimental observa- 
tions of the current may give important insight into var- 
ious fundamental properties of the system, as e.g. elec- 
tron correlations and the spin-charge separation. On the 
other hand, by controlling the rate of change of the ex- 
ternal magnetic flux, we demonstrate how an oscillating 
current of desired amplitude and frequency can be in- 
duced. Moreover, its time-average is non-zero and con- 
tains a dc component. The significant advantage of the 
method based on the magnetic flux variation is the 'non- 
invasive' manipulation performed outside the ring, with- 
out coupling to external leads. Recent progress in the 
highly controlled fabrication of quantum ring structures 
makes the verification of our findings quite realistic in 
the nearest future. We should note that basic limitations 
concerning mesoscopic systems (time shorter than the re- 
laxation time and the flux of order of the flux quantum) 
are not so restrictive for experiments performed in the 
optical lattice setup. In this context, our results pose 
some important implications for the design of quantum 
motors discussed in Ref. [lfj. For U/J < 1 or U/J > 1 
significant currents can be generated neither by impulses 
with <fr(tf) ~ 4>(U) nor by a magnetic flux that has small 
time-averaged value <f>[t) <C <f>a- Consequently, in the 
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